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ABSTRACT 

With its increased sensitivity and resolution, the Planck satellite is expected to improve the 
measurement of most cosmological parameters by several factors with respect to current 
WMAP results. The actual performance however, may depend upon various aspects of the 
data analysis. In this paper we analyse the impact of specifics of the data analysis on the ac- 
tual final results. We also explore the synergies in combining Planck results with future galaxy 
surveys. We find that Planck will improve constraints on most cosmological parameters by a 
factor 3-4 and on the tensor-to-scalar ratio r by a factor 9. Also inflationary parameters, like 
r, n s and n run , are practically not degenerate any longer. The tensor spectral index, however, 
is little constrained. A combination of the 70 to 143 GHz channels will contain about 90% 
of all possible information, with 143 GHz polarisation information carrying about half of the 
constraining power on r. Also, the error on r degrades by a factor 2 if no B modes are con- 
sidered in the analysis. High-/ temperature information is essential for determination of n s 
and Qb> while improving noise properties increase the /-range where Planck would be cosmic 
variance limited in polarisation, implying a significant improvement on the determination of 
r, t and A s . However, a sub-percent difference in the FWHM used in the data analysis with 
respect to the one in the map will result in a bias for several parameters. Finally, Planck will 
greatly help future missions like LSST and CIP reach their potentials by providing tight con- 
straints on parameters like n s and n run . Considering Planck together with these probes will 
help in breaking degeneracies between fix and fl\ or fi^m an d fv, resulting in improvements 
of several factors in the error associated to these parameters. 
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1 INTRODUCTION 

Observations of the cosmic microwave background (CMB) have 
driven an incredible improvement in our understanding of the Uni- 
verse in the last few decades. The discovery of CMB anisotropies 
by the COBE satellite triggered the planning of new space missions 
targeted to the study of CMB anisotropies: WMAP and Planck. 
The WMAP satellite has already delivered results and is still fly- 
ing, while the Planck satellite is scheduled to fly next year. It is 
therefore appropriate at this time to use the experience acquired 
with the WMAP data analysis in order to gain information on what 
can be expected in terms of parameter determination from Planck. 

Given the improved technical performances, Planck carries 
great expectations and several authors have proposed a plethora 
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lLa Vacca & Colombo! 2008) . In general, these estimates aim at 
an extended set of parameters beyond the minimal cosmological 
model, but rely on a number of simplifying assumptions. Actual 
constraints will depend sensitively on the ability of Planck to clean 
foregrounds from individual frequency channels, maximising the 
cosmological information available. 

In this paper, we adopt a different perspective focusing on 
quite simple cosmological models and then exploring the required 
Planck performance needed for useful constraints. We explore the 
effects of foreground cleaning restricting the frequency channels 
and angular scales available for cosmology, and the effects of beam 
degradation. These have differing effects on individual parame- 
ter constraints and weakened experimental performance can en- 
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hance degeneracies between parameter s. These issues have been 
partly addressed in the Planck bluebook {The Planck Collaboration 
2006). This paper is meant to provide further and more detailed in- 
formation than that contained in the bluebook, especially in light 
of what has been learnt from WMAP. At the same time, we hope 
that our results are more transparent than those of full experimen- 
tal simulations, which fully include experimental imperfections but 
are restricted in the range of parameter space explored. 

Specifically, we will first consider a minimal set of parame- 
ters and perform a Monte Carlo Markov Chains (MCMC) estimate 
of the accuracy Planck will have in determining them, also show- 
ing under which conditions these results are attained. Then we will 
consider a more extended set of parameters including the possi- 
ble presence of gravitational waves, curvature, neutrino-related pa- 
rameters, quintessence and we will discuss how the previous re- 
sults are modified in the presence of this extended set. While an 
appropriate treatment of Planck paramete r estimation should also 
contain a general approach to reionisation jMortonson & Hul2008l : 
IColombo & Pierpaolill2003) . in the following we will simply con- 
sider a sharp reionisation pro cess parametrised only by th e total 
optical depth r. As shown in IColombo & Pierp aoli (2008), such 
assumption should not impact the size of errorbars for most pa- 
rameters, which is the main interest of this paper. Finally, we will 
discuss the role of Planck as support for the success of galaxy sur- 
veys as cosmological probes. Galaxy surveys are able to constrain a 
subset of the cosmological parameter set, typically relying on CMB 
data to provide information on the remaining parameters and to sig- 
nificantly reduce degeneracies. Currently, no space CMB mission 
is planned after Planck, and therefore the combination of Planck 
with external datasets such as galaxy surveys will be a fundamen- 
tal tool for estimating the cosmological parameters in the coming 
years. 

The simplified procedure adopted here will by no means be 
the one adopted during the actual data analysis. For example, the 
effect of the precise sky cut appli ed on power spectrum estima- 
tion and likelihood evalu ation (e.g.. lHivon et alj|2002l ; lLewis et al.l 
l2002l ; lBrown et a l. 2005), th e effect of foreground residual o n small 
and large scale (e.g-. ISerra et alj|2008l ; iBetoule et alj 2008) or the 
effect of lensing on th e polarisation power spectrum are all not 
investigated here (e.g., lLewisll2005l ISmith e t al. 2006). All these 
extra complications will however tend to degrade the level of per- 
formance for the Planck satellite. The results presented here should 
therefore be considered as a target or an optimistic limit. 

This paper is organised as follows: in sectionf2]we will present 
the method and apply it to a configuration mimicking the WMAP 5 
yrs data release; in section[3]we will discuss the performances ex- 
pected for Planck for different frequency combinations and sets of 
parameters; section|4]is dedicated to investigating the role of Planck 
as support for other missions and finally section [5] is dedicated to 
the conclusions. 

2 METHOD AND WMAP-5YR CASE 

In this paper we will perform a MCMC analysis of simplified 
mock data to estimate Planck perfor mances in parameter es tima- 
tion, using the CosmoMC packagtQ l iLewis & Bridle! |2002|) The 
availability of fast codes for evalua tion of CMB power spectra, 
e.g. PICCq dFendt & WandertH2007r) . allows good convergence of 

1 http://cosmologist.info/cosmomc/ 
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the chains to be obtained in ~ lh on modem office worksta- 
tions, thus providing a valid alternative to a Fisher matrix approach, 
while avoiding so me of the Fisher matrix shortcomings (see, e.g., 
IPerotto et alj200rj , for a more in depth comparison of MCMC and 
Fisher matrix approaches). 

Specifically, we will consider idealistic simulations of CMB 
data, assuming white isotropic noise and Gaussian beams. Real data 
analysis will have to deal with complications like anisotropic and 
correlated noise, beam systematics, calibration effects. A complete 
treatment of these effects is beyond the scope of this paper and to 
some extent would require access to the actual measurements and 
in-flight data. In addition, we do not consider residual foreground 
contribution in the mock data. Rather we assume that some of the 
frequencies will be used to completely clean the remaining chan- 
nels, which will be the only ones used for data analysis. In pra ctice, 
this is the strategy adopted up to now (e.g. lDunklev et al . 2008). We 
will then show how our results depend on this choice of channels. 

For full sky and noiseless data, the exact likelihood of a cos- 
mological model, defined by a set of CMB angular power spectra 
C/ h , giv en a simulated dataset, Ci, is given by the Wishart dis- 
tribution iPercival & B rown 20061: lHamimeche & Lewisll2008h . It 
reads 

log£(C £ th |C f ) = (1) 

-i^(2^ + l){Tr[c f (C/ 11 )- 1 ] -lnlCMC/'T 1 ! - n) 

i 

where n is the number of different modes (T, E and B) consid- 
ered. With this choice of nor malisation log£(C/ h |C^) = for 
C/ h = Ci . In the presence of white noise and for an instrument 
with Gaussian beams, the above expression holds provided we re- 
place Ce with Ci +Af{Bi, where JVe are the white noise spec- 
tra and Bp is the spherical harmonic transform of the instrument's 
beam. In addition, when combining different frequencies, we use 
an inverse noise weighting. 

In the presence of sky cuts or non-uniform noise, equation Q] 
is no longer exact. In this work, we approximate the effect of sky 
cuts by multiplying the r.h.s. in equationQ]by a factor / s 2 ky , where 
/ak y is the fraction of sky actually observed. This ad-hoc correction 
crudely accounts for the loss of statistics and the induced correla- 
tions between different multipoles, and is already implemented in 
CosmoMC. For Planck, we assume / s k y = 0.80. This approxima- 
tion is not well suited to describe the behaviour of the likelihood 
function at low multipoles, which are critical for the determination 
of cosmological parameters like the optical depth to reionisation, 
r, and the tensor-to-scalar ratio, r. Features like non-uniform sky 
coverage alter the shape of the likelihood function, and introduce 
correlations between multipoles, in particular at £'s corresponding 
to scales similar to the dimensions of the cuts. A discussion of how 
this affect our results for r is presented in section [T2l 

As a first application, we consider the WMAP 5 yrs (WMAP5) 
case, introducing in the analysis the frequencies that were actually 
used in the real data analysis at different scales. In particular, we 
adopted V and W bands in temperature and Ka,Q and V bands for 
large scale polarisation reducing to Q and V only for small scale 
polarisation. For WMAP we considered 75% sky coverage, accord- 
ing to the mask applied by the WMAP team to the actual observed 
maps, and an average noise value for each frequency considered 
computed from such maps. We do not consider the effect of lens- 
ing or marginalisation over the amplitude of the Sunayev-Zeldovich 
template. 

The results are summarised in table Q] where we report the 
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Table 1. Comparison of uncertainties from actual WMAP 5-year data 
analysis and MCMC estimates. Table shows the percentage difference 
1 — ff MCMc/ ff WMAP f° r the basic ACDM + tensor model. The third 
column shows the MCMC estimates of the standard deviation in unit of the 
expected Planck accuracy, assuming cleaning of the 70, 100 and 143 GHz 
channels (see table[3J. 
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Table 2. Planck specification as reported in the Planck bluebook. Listed sen- 
sitivities are goal sensitivity assuming 14 months integration, for a square 
pixel with a side equal to the angular resolution of the corresponding chan- 
nel. 



percentage difference in the parameters' sta ndard deviation with 
respect to the actual WMAP5 performances iDunkle v et alj[2008h 
for a ACDM cosmology with tensors. For comparison, we also re- 
port the expected performance of the Planck satellite, considering 
an ideal combination of the 70, 100 and 143 GHz channels (see 
below for more details). We conclude that this simplified approach 
in the case of WMAP5 leads to underestimate of parameters uncer- 
tainties by 25 — 30%. Although this result does not fully validate 
this procedure for Planck, it is encouraging that its performances on 
simulated WMAP5 data are so close to the actual measured value. 

FigureQjprovides a visual comparison of our analysis applied 
to simulated WMAP5 and Planck data. Planck is expected to im- 
prove on this result by typically a factor 2.5 to 4, except in the case 
of r, whose upper limit for a fiducial model with r — shrinks by 
almost an order of magnitude. It is also evident that Planck will re- 
move or strongly reduce several degeneracies which affect WMAP 
data, in particular those involving the optical depth to reionisation, 
r. Notice that in the WMAP5 simulations we did not include the 
contribution of B-modes, while for Planck we also included sim- 
ulated BB data. If systematics and/or foregrounds prevent Planck 
from measuring the cosmological BB spectrum, the expected im- 
provement on the tensor-to-scalar ratio upper limit would be by a 
factor ~ 3 (see section [3~2l for a detailed discussion). 



3 PLANCK-ONLY CONSTRAINTS 

Planck has a noise level which is a factor 2-10 lower than WMAP5, 
a wider frequency coverage, and a finer angular resolution. Cosmo- 
logical constraints for Planck will therefore come from a different 
set of scales and frequency ranges than in the case of WMAP. 

In the following, we will adopt the simplified procedure out- 
lined above to assess Planck performances under different data 



analysis circumstances. We will use the noise specifications derived 
from the bluebook, also reported in table [2] 

Initially we will consider the following fiducial ACDM min- 
imal model, based on WMAP5 results: flat, cj b = 0.02273, w c = 
0.1099, r = 0.087, n s = 0.963, H = 71.9 Km/s/Mpc and no 
tensor perturbations. The amplitude of the scalar fluctuation spec- 
trum, A a , derives from the requirement that as ~ 0.796. Following 
CosmoMC convention, we define here the tensor-to-scalar ratio r 
as the ratio of the initial curvature and gravitational waves power 
spectra at a reference scale of k plv = 0.05Mpc _1 . In the rest of 
the paper, we define .4 S = Log(10 10 yl s ). In order to explore the 
parameter space in a more efficient way, it is useful to chose as 
a primary parameter the angle subtended by the sound horizon at 
recombination, 9. Our basic parameter set is then: 

{uj b ,uj c ,9,r,n B ,A B ,r} . (2) 

Notice that we include r in the analysis, even if the fiducial model 
has no tensors (see section [Jl2l for a discussion of models with non- 
zero r). With this choice of parametrisation, Ho and as are then 
derived parameters. 

Plots for power spectra in this minimal model are compared 
to the noise level and cosmic variance at different frequencies in 
figure [2] In temperature, all Planck channels will be "cosmic vari- 
ance limited" (i.e. the noise equals the cosmic variance) in each 
frequency channel up to angular scales I ~ 800; while the combi- 
nation of all channels up to 217 GHz is cosmic variant limited up 
to £ ~ 1000. 

The channels for which the signal-to-noise ratio (S/N) is unity 
at the smallest scale (of £ ~ 1600) are 143 and 217 GHz, the former 
having better signal to noise at larger scales and the latter at smaller 
scales. Combining all channels, the signal-to-noise ratio becomes 
unity at I ~ 1800. These figures partly depend on the values of the 
reference cosmological parameters, in particular A s . WMAP5 data 
constrain A s at < 10% level; the effect of this uncertainty on the 
above discussion is negligible. 

As for polarisation, Planck will be cosmic variance limited up 
to i ~ 8 - 12 (cfr with t = 4 - 5 for WMAP5), depending on the 
cosmological model considered and, in particular, on the details of 
the reionisation history. The smallest scale with S/N above unity is 
I ~ 800 for a combination of channels up to 143 or 217 GHz. 

In the following, we will explore how parameter constraints 
would vary in the case that only a subset of frequency channels 
could be used for extracting cosmological information. 

3.1 Minimal model: results on parameters 

A great advantage of multi-frequency experiments is to be able 
to use frequency information to subtract foreground contributions 
from the maps and combine the information of different channels 
to derive cosmological parameters. In the case of WMAP5, for in- 
stance, a great deal of leverage in parameter estimation has come 
from the ability to use the Ka band in polarisation on large scales. 

Similarly, Planck's performance will depend upon the ability 
to clean foregrounds from the largest number of channels. Planck 
will be able to capitalise on the existing WMAP data for its lowest 
frequency channels, but being an experiment with higher sensitiv- 
ity it may need to use its own channels to perform the job fully. 
The most crucial concerns, of course, arise regarding polarisation, 
which is currently poorly measured. As there is consensus on the 
fact that polarised foregrounds are mi nimal around 70 GHz (e.g. 
ISamtleben et"afl2 007; Gold et al. 2008), we will consider as a min- 
imal assumption that at least the 70 GHz channel will be cleaned. 



© 0000 RAS, MNRAS 000, 000-000 



4 L.P.L. Colombo et al. 








0.2 0.4 



Figure 1. Marginalised distributions and joint 2D confidence regions for the base fiducial model analysed assuming a combination of 70, 100 and 143 GHz 
Planck channels (blue/cyan lines) or WMAP 5yr specifications (red/orange lines). The 2D plots show the 1- and 2-<r regions; crosses mark the input values 
for the parameters. 
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Table 3. Error estimates for different combinations of Planck channels for the reference ACDM + tensors model. For all parameters we report the 
standard deviation for the marginalised distribution, except for r for which we report the upper 95% confidence limit. Second column: estimates 
assuming cleaning of 70, 100 and 143 GHz channels. Values shown are actual errors. Columns 3-8 shows estimates for different combinations of 
channels normalised to the values of column 2, when a frequency is followed by a T we only consider temperature data for that channel. 
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WMAP5 was not able to use the 95 GHz channel in polarisation 
due to potential dust contamination, however, thanks to the greater 
frequency coverage of Planck, we will assume that the 100 GHz 
channel will also be available to extract cosmological information. 
Starting from this minimal configuration (70 + 100 GHz), we will 
then explore the effect of adding information from higher and lower 
frequencies, in temperature only or both in temperature and polar- 
isation. The results are summarised in table[3] where we quote the 
estimated errors on different parameters with respect to the case for 
the combination of channels 70 + 100 + 143 GHz. This frequency 
combination is, in fact, already providing most of the information, 
so that adding all LFI frequencies and the 217 GHz channels would 
improve the results by only 10% in all parameters. The higher fre- 
quency channels would add even less, due to their higher noise 
level. Not being able to use the 143GHz channel in polarisation 
increases errors by about 15 — 25% (n s being the least impacted), 
except for r in which case the upper limit roughly doubles. This 
problem is mitigated if all LFI channels are present. Not being able 
to use the 143 GHz channel at all, however, would degrade the cos- 
mological parameter determination by 20-50% on most parameter 
and a factor 2 on r. Also in this case, using all LFI gives ~ 10% 
improvement on r and A s while other parameters are affected at 
the percent level. 

This simplified analysis summarises the information con- 
tained at all scales. It is indeed the case that certain foregrounds 
are a concern only for either temperature or polarization, or only 
on some spatial scales. One such example is residual point sources, 
which may spoil the use of small scale temperature fluctuations. It 
is therefore interesting to ask how much of the information comes 
specifically from temperature measurements and which are the rel- 
evant scales in determining each parameter. In order to address this 
question, we considered the 70 - 217 GHz channels combination 
and evaluate how much each error estimate decreases as a function 
of the maximum multipole £ max considered in the analysis. Re- 
sults are compared to the case of an ideal experiment with cosmic 
variance limited (CVL) temperature measurements over the whole 
range of multipoles considered, and polarization sensitivity equal 
to that of the selected Planck channels. The output of this analysis 
is reported in figure[3] 

In principle, for the minimal cosmological model considered 
here, Planck will recover all the relevant information on parameters 
encoded in the CMB power spectra up to scales I ~ 1500, after 



which differences with a CVL experiment start to show. The pa- 
rameters whose determination is most improved when temperature 
information on scales above £ — 800 is included are n s , uib and 8 
(by about 30-40%) while r, A 3 and r are basically unaffected. 

While Planck will perform, in temperature, as well as a CVL 
experiment up to I ~ 1500, a greater room for improvement is 
left for polarization measurements. This improvement may come 
either if Planck keeps functioning after the 14 months required to 
complete 2 full sky surveys, or by other ground-based or balloon- 
borne experiments, like the currently planned SPIDER dCrill et all 
2008). SPIDER will cover ~ 50% of the sky, with a polarization 
sensitivity higher than Planck, but with a lower angular resolution. 
It is then interesting to consider how the combination of Planck 
temperature measurements with better polarization data will affect 
cosmological parameter estimation. To asses this point we consider 
an ideal experiment with temperature information equal to the com- 
bination of the 70 — 217 GHz Planck channels and CVL polariza- 
tion measurements up to a multipole lev ranging between £ — 10 
and £ — 800. Above £cv polarization sensitivity is equal to that 
expected from Planck. Let us also point out that CVL determina- 
tion of i5-mode polarization above I of a few hundreds will likely 
require the next generation of spa ce CMB mission, like the recently 
proposed EPIC (Boc k et al . 2008). However, if such a mission does 
not have high angular resolution, we do not expect a significant im- 
provement on Planck temperature data, which are effectively CVL 
for I < 1000. 

The right panel of figure[3]shows the improvement on param- 
eter constraints over Planck 70 — 217 GHz, as a function of £cv- 
As expected, a CVL determination of even the first 10 polarization 
multipoles would significantly improve the constrain on r. Results 
shown refer to a fiducial value r = 0.05, for which lev = 10 
would improve Planck constraints by a factor ~ 2; the improve- 
ment would be even more relevant for lower tensor-to-scalar ra- 
tios. While the large angle B-mode spectrum is expected to be 
at most comparable to Planck noise, the i5-mode spectrum will 
be measured with essentially a CVL accuracy up to I — 8 — 9, 
so that lev = 10 will provide only a modest, ~ 10%, improve- 
ment in the determination of r. This improvement will increase to 
~ 40% if CVL data extend up to £ ~ 20; above this value the 
polarization C\ are mostly independent of r for the class of sharp 
reionisation history considered here, and further improvements on 
sensitivity do not lead to significant improvements on determina- 
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Figure 3. Le/f.Sensitivity to the different parameters as a function of the maximum multipole considered in the analysis. Plot shows error estimates for the 
various parameters as a function of ^ max considered in the analysis, normalised to the error for ^ max = 800. Heavy lines are for the combination of the 
70, 100, 143 and 217GHz Planck channels, thin lines are for an ideal experiment with cosmic variance temperature measurements and polarization sensitivity 
equal to the Planck configuration considered here. It is clear that up to <? max ~ 1500 a CVL experiment would not offer a significant advantage over Planck. 
Right : Impact of higher sensitivity polarization data on parameter constraints. Plot shows parameter estimates for an ideal experiment with CVL polarization 
data up to lev and sensitivity equal to Planck 70 — 217 GHz channels for I ^ lev- Error estimates are normalised to Planck estimates, corresponding to 
<?CV = 2. 



tion of the optical depth. However, the multipole range 20 — 50 
contains information on reionisation for non-minimal reionisation 
histories, and improving determination of this part of the spectra 
can significantly increase constraints on these models. Due to the 
r — A a degeneracy, the normalisation of the scalar power spectrum 
improves in a similar manner, which in turn affects the determina- 
tion of as. A second set of parameters, including co c , n s and Ho, 
shows a ~ 15% improvement for lev = 100 — 200, increasing 
to ~ 25% for lev = 800. A third group of parameters, 9 and u/b, 
shows only a modest improvement for CVL data extending to the 
first 500 multipoles, and a ~ 20% improvement for lev = 800. 

In what follows, whenever not specified, we assume the fre- 
quency combination 70 + 100 + 143 GHz and include in the anal- 
ysis multipoles up to £ max = 2500, even if for Planck multipoles 
I > 2000 only have a (few) percent level impact on parameter de- 
termination. 



3.2 The tensor-to-scalar ratio determination 

As it appears from table Q] one parameter whose determination 
Planck is expected to greatly improve is the tensor-to-scalar ra- 
tio r. For a fiducial model with negligible tensor contribution, r < 
0.01 corresponding to s mall-field inflation m odels (for a review of 
inflationary physics, see lLvth &~R iotto 1999), using the 70 — 143 
GHz Planck channels we expect an upper limit r < 0.03(95%c.L), 
tightening WMAP5 constraints by about 1 order of magnitude. 
Without the polarization information of the 143 GHz channel, this 
constraint degrades by more than a factor 2, while adding the 217 
GHz channel would improve this figure by ~ 15%. If it proves 
possible to use all the LFI channels for cosmological parameter es- 
timation, the upper limit would improve by ~ 25% if the 143 GHz 
polarization is not included in the analysis, or by ~ 10% otherwise. 

For a non-negligible, i.e. r > 0.05, contribution by gravita- 
tional waves, Planck will be able to put a lower limit on r with a 
confidence level > 95% (see figureQ. It is then interesting to ex- 
plore to what extent results depend on the frequencies considered. 
Results for a fiducial cosmological model with r = 0.05 are given 
in table [4] A combination of the 70-143 GHz channels allows r 



to be constrained with an accuracy of ay ~ 0.023. The uncertainty 
on r increases by ~ 50% without the polarization information from 
the 143 GHz channel, however, since the marginalised distribution 
is markedly non-Gaussian even in this case it is possible to claim a 
2cr detection of r > 0. Adding the 30 and 40 GHz channels would 
improve constraints on r by 10 — 15%. Notice that these results 
were obtained by fixing the tensor spectral index to tit = (see 
below for further discussion). 

The BB spectrum is a unique signature of tensor perturba- 
tion, however it is also a weaker signal than the EE and TE spec- 
trum; it has not yet been observed, and is potentially more affected 
by foreground re siduals. Simulation s of polarised foregrounds re- 
moval strategies (Betoule et alj2008h suggest that for Planck resid- 
ual foreground contamination would increase the uncertainty on r 
by ~ 30%, for a fiducial value r — 0.10. Thus, we expect fore- 
ground cleaning to have an impact on the uncertainty on r which is 
(at least) comparable to that of the instrumental noise. However, a 
detailed knowledge of polarised foregrounds is still lacking, and a 
correct assessment of Planck capabilities of measuring the B-mode 
Ci will require the actual data. 

For these reasons, it is sensible to ask how much of the con- 
straining power resides in the BB measurement and to what extent 
the constraint on r would be weakened if instead we were only able 
to measure TT, EE and TE. The results are summarised in figure [4] 
for r = 0, 0.05, 0.1, 0.15, 0.2, and a fiducial set of parameters as 
above. In the case of a high value of r that Planck might possibly 
detect, not being able to use the BB spectrum would imply a dou- 
bling of the error bars, and a 3<r detection would be possible only 
for va lues r > 0.15 — 0.20, in agreement with iKnox & Turner! 
dl994h . 

It is then clear that a significant portion of Planck's capabil- 
ities of detecting tensor modes comes from the \ow-l part of the 
BB spectrum. When dealing with real data, contamination by the 
Galaxy, as well as component separation and point sources sub- 
traction, will require masking of parts of the sky. In this case, the 
correct shape of the likelihood function is no longer given by equa- 
tion Q] and for arbitrary cuts it is not possible to write an analyt- 
ical expression for the likelihood in harmonic space. As discussed 
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above, we assumed here that it is possible to account for sky cuts 
simply by multiplying the Wishart distribution by a factor /JL.. To 
test that this approximation does not significantly affect our deter- 
mination of r, we compare it with the exact likelihood function 
evaluated in pixel space (we refer the interested reader to, e.g., 
Tegm ark & de Oliveira-Costall200lL for the relevant expressions). 
Evaluating the full Planck likelihood in pixel space is not numeri- 
cally feasible, since it would require the inversion of a matrix with 
a side of O(10 7 ) elements. However, we focus here on the determi- 
nation of r which depends mainly on the I < 30 multipoles, so that 
we ca n work with low resolution maps. Using the HEALPix pack- 
ag<0 ( Gorsk i et alj20ol) , we generate temperature and polarization 
maps at a resolution iV s idc == 16, corresponding to a pixel size of 
~ 3.5, including contribution by multipoles up to I = 32. We add 
a white noise corresponding to that expected for the 70 — 140 GHz 
channels; due to the low resolution of the maps, the beam finite res- 
olution is not relevant. From each map we remove the pixels with 
centers within 10 of the equator. Using the remaining pixels, cor- 
responding / a ky = 0.8125, we evaluate the exact likelihood for r, 
with other parameters fixed to the fiducial value. We compare this 
likelihood to the one evaluated with equation [TJ using as fiducial 
spectrum the full-sky spectrum of the CMB + noise map. Both in 
pixel and harmonic space, we include in the likelihood evaluation 
TT, TE, EE and BB multipoles up to £ = 32. We repeat this test 
for two fiducial values r — and r — 0.10, for each value we 
perform 50 sky realization. For r = we compare the 95% up- 
per limits of the pixel and harmonic based likelihood, and find that 
in average the approximate likelihood we use throughout this work 
overestimates the correct result by ~ 4%. For r = 0.10, we instead 
compare the variance of the two sets of distributions; using the ap- 
proximate likelihood overestimates the correct results by ~ 11%. 
As expected, when considering single realizations, the estimate of 
r obtained with the rescaled Wishart distribution are different from 
that obtained in pixel space, as the former likelihood peaks at the 
full-sky power spectrum, while the latter correctly takes into ac- 
count only the observed portion of the sky. On average this dis- 
crepance is ~ 0.35oy. Elsewhere in this work we build the mock 
data using the theoretical power spectrum, instead of considering 
realizations, therefore there is no bias in the corresponding results. 
Since r is not significantly degenerate with the remaining param- 
eters, with the exception of ut, we expect that these conclusions 
hold also when fitting for the remaining cosmological parameters, 
and that results of this paper do not significantly depend on the sky 
coverage of actual data. 



3.3 Beam degradation 

Planck results will come from a combination of high sensitivity and 
fine angular resolution. While a comparison of Planck with a CVL 
experiment gives an assessment of the relevance of high sensitiv- 
ity, a crude test of the impact of angular resolution can be made 
by comparing Planck's expected performance with that of an ideal 
experiment with ad-hoc increased beam FWHM at all frequencies 
by either 10% or 30%. 

Table|5]reports the parameter uncertainties for such an experi- 
ment, relative to the Planck expected performances quoted in table 
[3] An increase of 10% in FWHM worsen the constraints on most 
parameters by ~ 5% — 10%, while degrading the beams by 30% 



3 http://healpix.jpl.nasa.gov/ 
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Table 4. As tablef3] but assuming a non-zero tensor amplitude, r = 0.05. 
In this case for r we quote the standard deviation as done with all the other 
parameters. 
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Table 5. Error estimates for an experiment with sensitivity equal to Planck 
70-143GHz channels and larger beams. Values are quoted relative to Planck 
estimates reported in table [3] 



increases error estimates by ~ 10% — 20%, for our basic config- 
uration. As shown in figure [5] some parameters, most noticeably 
n s , ujb and 9, gain constraining power from the measurements of 
high-^ TT power spectrum. 

Real beams can be affected by a number of systematic effects, 
e.g. deviations from Gaussianity or asymmetries. A full assessment 
of these effects will require in-flight beam calibration, and is out- 
side of the scope of this paper. Here instead, we perform a simple 
test of beam impact on parameter estimation by assuming that the 
Gaussian FWHM used in the analysis is different (either bigger or 
smaller) by the actual FWHM by 0.05% or 0.20%, for all frequen- 
cies considered. The major impact of a beam mismatch is a bias on 
most cosmological parameters, as shown in figure|6] Even an error 
of 0.05% in the assumed beams leads to a bias of ~ 0.5a on most 
parameters, increasing to more the ~ 2a for a misestimate of the 
FWHM by 0.20%. For Gaussian beams, the beam transfer func- 
tion depends on a single I. Thus, if the beam mismatch becomes 
relevant only at scales I > 1800 — 2000, it is expected to affect 
parameter estimation at the few percent level (see figure [3). Real 
systematics, however, will also introduce correlations between dif- 
ferent multipoles so in general it will not be possible to reduce these 
effects by simply excluding this range of £'s. 



3.4 Away from the minimal model 

After having explored Planck's performance for the minimal cos- 
mological model, we now broaden the parameter space and investi- 
gate performance for a series of less minimal but still quite general 
cosmological models. These will include a running spectral index, 
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TT+TE+EE+BB blue 



i l I 
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Figure 4. Impact of 73-modes detection on the determination of the tensor-to-scalar ratio, r, for different fiducial values of r. For r > bars show the 
expected lcr error on the tensor-to-scalar ratio including (blue) or excluding (gold) B-modes in the analysis. For a fiducial r = bars show the upper 95% 
c.l. instead. If B-modes information is not available, error on r increase by a factor ~ 2. Notice that for r < 0.10, the corresponding marginalised distribution 
is markedly non-Gaussian and quoting the standard deviation does not properly characterise the la confidence interval. Results shown are for the 70 — 143 
GHz channels combination. 
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Figure 5. Impact of small angular scales on parameter determination. Histograms shows the uncertainties on uj-^, n B and 9 for nominal beams (left bar), beam 
enlarged by 10% (middle) or by 30% (right). For each parameter, the left panel refers to the 70 — 143GHz channels, the right panel is for 70 — 217 GHz 
channels. 



neutrino-related parameters, curvature and quintessence. For these 
models, we will explore frequency dependence and degeneracies. 

3.4.1 Dark Energy and Cosmic Curvature 

In most cosmologies, the size of the sound horizon at last scatter- 
ing provides a standard ruler at a redshift of zls — 1090. Indeed, 
most information on the geometry of the Universe present in the 
CMB power spectra, can be summarised in terms of the redshift 
of photon decoupling, z„, and of tw o distance ratios, the acoustic 
scale I a and the shift parameter R dElgar0v & Multamakil l2007l ; 
Mukh eriee et alj 2008; Ko matsu et al.l 120081) . In principle this al- 
lows for a substantial compression of information without a signif- 
icant loss of accuracy, especially when combining CMB data with 
external measurements. This approach however cannot be adopted 
when simultaneously fitting for the entire set of cosmological pa- 
rameters we are considering. 

The geometrical degeneracy limits how the accuracy on the 
spectra translates in accuracy on the cosmological parameters 
defining the geometry of the Universe or dark energy properties and 
current constraints on these parameters rely heavily on external in- 
formation. We wish to investigate here to what extent this will be 
true for Planck as well. In particular, we will consider cosmic cur- 
vature, Qk, and a dark energy component with a constant equation 
of state, w. 

Constraining the geometry of the Universe to high precision 
provides a weak, but important probe of inflation. Generic inflation 
models tend to produce a Universe completely flat except for super- 
horizon curvature fluctuations at the level of Q,k ~ 10~ . Although 
this high level of precision is likely beyond currently envisaged ex- 
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Table 6. As table [5] but for a ACDM + tensor + curvature model. 



periments a positive detection of curvature at a level higher than 
this would be a very interesting result. This would not necessarily 
invalidate the inflationary p aradigm, since open inflation models 
dLinde & Mezhl umiarJll995h can produce Ok S; 10 3 . Measure- 
ments at this level would move towards constraining such models. 

Table [6] reports results when the parameter space is expanded 
to include Q,k- Our default Planck configuration will constrain cur- 
vature with an accuracy of an K ~ 0.023, compared to an accuracy 
of an K ~ 0.050 for the WMAP5-like case (see figure[7]l. Allowing 
for curvature degrees of freedom in the analysis increases the error 
on Hq by more than an order of magnitude and almost doubles 
the error on as ■ As expected, other parameters are not significantly 
affected. 

We now turn to discussing estimation of dark energy prop- 
erties. The simplest way to account for dark energy is to ad-hoc 
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Figure 6. Effect of an incorrect beam characterisation on parameter estimation for Gaussian beams. Figure shows the bias in the recovered parameter resulting 
by underestimating (red lines) or overestimating (green lines) the actual beam FWHM by 0.05% (dashed lines) or 0.20% (solid). 




Figure 7. Joint 1— and 2 — a confidence regions in the S1a — (left) and in the — fip; (right) planes. Plots compare results for the 70 — 143GHz 
channels (blue/cyan lines) with those for simulated WMAP5 data (red/orange). Crosses mark the input value of the parameters. Planck will significantly reduce 
the degeneracies present in WMAP5 data. 



include a cosmological constant term into Einstein's equations. 
A more physically moti vated alternative is provided by a self- 
interacting scalar field ilRatra & Peebles! 1 19881 : IWetterichl 1 19881 ; 
ICaldwell et afll 19981 : iBrax & Martij|l999h . A common prediction 
of these models is that dark energy can be described as a perfect 
fluid with a time varying equation of state, w(a). In general, the 
effects of dark energy on CMB spectra are well described by mod- 
elling dark energy as a perfect fluid with constant equation of state 



J da Sl w (a)w(a) 
J da fin, (a) 



(3) 



This approximation is not valid for models with direct dark energy 
- dark matter interactio n (e.g., lAmendola|[2000l : iHuev & W andelt 
20061; iMainini et ai]2 005 ) and may lead to si gnificant biases in th e 
determination of cosmological parameters iVergani et al.l l2008h . 
We will not consider the latter class of models here, and parametrise 
dark energy with a constant equation of state, w. We also assume a 
flat Universe. 

Error estimates for WCDM models are quoted in table[7] while 
figure[8]compares degeneracy regions for Planck and WMAP in the 
w — and w—ag planes. As with curvature, allowing for w ^ — 1 
significantly increases the errors on Ho and ag. In this case, how- 
ever, estimates are not dependent on the combination of channels 
used, and in addition Planck provides only a marginal improvement 
over WMAP5. Thus, at Planck sensitivities constraints are driven 
by the CMB internal degeneracies and the prior used in the analysis 
(e.g. on Ho), rather than by the data. 

In summary, Planck will significantly improve over WMAP 
5-year constraints on Q,k, but not on w. However, Planck will sig- 
nificantly improve over WMAP5 if we are considering models with 
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Table 7. As table[3] ACDM + tensor + w model. Errors on to, <rg and Hq 
do not depend on the channels used. 



both non-zero curvature an d w ^ —1, or models with more than 
one dark energy parameter JXia et al.l2007h . In both cases, external 
data set will still be fundamental (also see Mukheriee et al. 2008). 



3.4.2 Neutrinos and relativistic degrees of freedom 

The detection of oscillations of solar and atmospheric neutrinos 
has confirmed that v's are massive particles. However, for cur- 
rentl y allowed values of n eutrino's mass, 2~2 m v IS 1.3 eV (95% 
c.l.) komatsu et alj|2 008), neutrinos will still be relativistic at the 
epoch of last scattering, so that the effective matter-to-radiation 
ratio at zls is lower than the matter-to-radiation ratio today. In 
turn, this results in a shift of the position of the first acoustic peak 
due to the faster decay of the gravitational potential around re- 
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Figure 8. Joint 1— and 2 — a confidence regions in the to — S1a (left) and in the w — <rg (right) planes, for a combination of 70 — 143 GHz Planck channels 
(blue/cyan lines) or for WMAP5 (red/orange lines). Cross marks the input value of the parameters. Results are dominated by intrinsic CMB degeneracies and 
by the choice of priors (especially on Ho)- 



combination, and alter the shape of the acoustic peaks dPierpaolil 
l2003l : llchikawa et al J 120051 iLesgourgues & Pastor! l2006h . In addi- 
tion, neutrino mass slightly affect the expansion rate at late times. 

The three standard neutrinos may not be the only free stream- 
ing relativistic particle species present in the early Universe. Addi- 
tional relativistic degrees of freedom, usually parametrised in terms 
of effective number of neutrino species, 7Y e fr, would push the red- 
shift of matter-radiation equality to later epochs, which in turn af- 
fects the ratio of the heights of the first and third acoustic peaks 
of CMB spectra. Moreover, extra relativistic degrees of freedom 
change the expansion rate at early times, in particular affecting the 
time between (standard) neutrinos decoupling and the opening of 
the Deuterium bottleneck, which fixes the abundance of light ele- 
ments produced in the Big Bang Nucleosynthesis. Knowledge of 
the Helium abundance, Yn c , would then tighten constrain on N e g. 
WMAP5 data alone are not able to place an upper limit on N e s, 
even assuming Yho is known. With its improved characterisation 
of the acoustic peaks, Planck i s expected to significa ntly constrain 
both these parameters (also see Ichi kawa et alj |2008). 

Here, we discuss the determination of neutrino masses sep- 
arately from that of the effective number of relativistic species 
and Helium abundance. We show how these constraints depend on 
the frequencies used, and how errors on the non-neutrinos related 
parameters increase with respect to results for a ACDM model. 
The fiducial cosmological model does not include massive neutri- 
nos, /„ = 0, or additional massless species, N e g = 3.04, and 
y Hc = 0.24. 

When expanding the parameter space to include v's mass, we 
assume that the three mass eigenstates are completely degenerate 
and parametrise the neutrino contribution to the total dark matter 
energy density by /„ = fi„,o/fidm,o with Q d . m ,o = ^c,o + n„, . 
Results for the 70 — 143 GHz are reported in table [8] We find an 
upper limit /„ < 0.072 (95% c.l.), which translates into a con- 
straint on the sum of v 's masses ^ m„ < 0.77 eV, according to 
Y.m v = 94 eVf2„,o/i 2 = 94 eV/„Q d m,o/i 2 - Using only the 70 
and 100 GHz channels, this limit increases by ~ 20%, while if the 
217GHz channel is included the limit improves by ~ 5%3 Allow- 



4 Due to the stochastic nature of MCMC methods, the la upper limit on 
/„ is determined with a 5% accuracy. 
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Table 8. Estimated accuracy for a ACDM + tensor + massive neutrinos 
model. We assume 3 families of massive neutrinos with fully degenerate 
mass eigenstates, and parametrise the massive neutrino contribution to the 
total energy density by f v . For f v we quote the 95% upper confidence 
limit. In parenthesis we report the ratio between estimates for this model and 
estimates for the base ACDM +r model for the same channels combination. 
As usual, results of columns 3 and 4 are quoted in units of uncertainties for 
the 70 - 143 GHz channels. 

ing for a contribution by massive neutrinos more than doubles the 
uncertainty on Wdm and increases errors on Ho by ~ 3.5 times. 
Since erg is an integrated quantity with a significant dependence on 
/„,Ud m and Ho, these effects combines so that the corresponding 
uncertainty increase by a factor ~ 7 — 8. Errorbars on other pa- 
rameters are only moderately affected, increasing by ~ 5 — 10%, 
implying no significant degeneracies with /„. When fitting for neu- 
trino properties, we assume a constant equation of state dark energy 
with w = — 1. This assumption does not affect our results since the 
m v - w degeneracy is relevant only when combining CMB data 
with external measurements such as galaxy surveys or Lya forest 
data {kannestad 2005; Komatsu et al]|2008h . unles s we also allow 
for a direct dark energy - dark matter interaction dLa Vacca et al.l 
l2008h . 

In addition to the increase in uncertainties, the degeneracy be- 
tween /„ and ojdm leads to overesti mate the dark matter density by 
~ la, as shown in figure [9] (also see lPerotto et alj2006t) . 

Fitting for the relativistic degrees of freedom and Helium 
abundance significantly degrades the accuracy on most parameters, 
in particular on uidm, d and n s , as shown by the results of table|9] 
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Table 9. As table ACDM + tensor + jV e jf + Yhc- Assuming cleaning 
of the 217GHz channel would increase accuracy on most parameters by 
~ 10%. The third column reports the errors in units of the uncertainty for 
the base ACDM +r model, see table|5] 
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Table 10. Error forecasts for the ACDM + tensor + running model. 



In agreement with lHamann et all J2008h and llchikawa et all 12008). 
we find that Planck alone will be able to constrain N e g to within 
~ 10%, even when simultaneously determining Ykc- 



3.4.3 Running spectral index 

A power law is the simplest possible parametrisation for the power 
spectra of both scalar and tensor perturbations. However, several 
inflationary theories predicts small deviations from a pure power 
law, typically expressed as a logarithmic running of the spectral 
index, n run = dn s /dlog(k), at a suitable reference scale, fc p i v . 
Here we choose k p w — 0.05Mpc _1 . 

If only the 70 and 100 GHz Planck channels are available for 
cosmological parameter estimation, we find a moderate degener- 
acy between n s and n run , similar to the one observed in actual 
WMAP5 data. This degeneracy results in an increase of the uncer- 
tainty on n s by ~ 10 — 15% over what is obtained for our min- 
imal parameter set. Including the 143 GHz channels substantially 
improves the determination of the high multipoles, which allows 
Planck to completely break this degeneracy, as shown in figure [lOl 
Even with the inclusion of the 143 GHz channel, a minor degener- 
acy between n Iun and r remains, which partly propagates to As and 
r. This degeneracy however, only slightly affects the constraints on 
these parameters. The remaining parameters are not significantly 
affected by the inclusion of n Iun in the parameter set, regardless of 
the combination of frequencies considered. 



3.4.4 Tensor spectral index 

One of the most ambitious goals of Planck and other future CMB 
mission is the detection of primordial gravitational waves. In the 
minimal cosmological model considered here, with the tensor spec- 
tral index fixed to the input value «t = 0, r is not significantly 
degenerate with the other cosmological parameters. However, even 
if the tensor-to-scalar ratio is quite high r ~ 0.15 — 0.20, Planck 
will be able to constrain at most the first 10 — 15 multipoles of the 
B mode spectra. It is then interesting to determine whether Planck 
will be able to simultaneously constrain r and ut- 

We consider two fiducial models with r — 0.01 and 
r = 0.10, and tit determined according to the single field 
inflation consistency relation tit ~ —r/8. In both cases, we 
also consider a non-zero running of the scalar spectral index 
n run = —0.02. The parameter set characterising our model is thus 
A s ,riT,r}. Results of the analysis are sum- 
marised in table II II 

We find that for Planck the normalisation and spectral index 
are completely degenerate, as shown in figure QT] Allowing for 
tit completely disrupts the instrument's capability for measuring 
r, even for moderate values of r = 0.10, as Planck does not have 
sufficient leverage on the B-modes to simultaneously constrain 2 
tensor mode parameters. Errors on the remaining parameters are 
similar to those shown in table [10] implying that m is not signif- 
icantly degenerate with the other parameters, in particular with n s 
and n Iun . 



© 0000 RAS, MNRAS 000, 000-000 



12 L.P.L. Colombo et al. 




Figure 10. Left: 1- and 2-a confidence regions in the ra s — n run plane, using 70 — 100 Ghz channels (dashed lines) or 70 — 143 GHz (solid lines) channels. 
Using just the 70 and 100 GHz channels allows for a mild degeneracy between the two parameters, which is broken if it is possible to clean the 143 GHz 
channel. Right: Confidence regions in the t — Umn plane. The contribution of the 143 GHz frequency does not help in breaking the mild degeneracy between 
the two parameters. 




Figure 11. Degeneracy between the tensor-to-scalar ratio, r, and the tensor spectral index, rat, for fiducial values {r = 0.10, rat = —0.0125} (left) and 
{r = 0.01, rat = —0.00125} (right). Even for moderate values r ~ 0.10, Planck does not have sufficient leverage on the B-modes to simultaneously 
constrain r and nt . 



4 PLANCK AS SUPPORT FOR OTHER MISSIONS 



In the previous section, we have seen that Planck's performance 
for wider parameter spaces is often limited by degeneracies in- 
herent to the CMB data. Further improvement on parameter con- 
straints will come not from improved CMB data, but by com- 
bination with other experimental techniques. In this section, we 
will use the Fisher matrix formalism to illustrate constraints for 
proposed galaxy surveys in combination with WMAP or Planck. 
The Fisher matrix technique h as bee n discussed in detail by many 
authors, e.g. Eisenstei n et alj dl999l). and we direct the interested 
reader to Pritchard & Pierpaolil d2008h for further information on 
our Fisher matrix methodology. Throughout we will assume usage 
of the Planck's 70-143 channels for cosmological constraints and 
consider WMAP5 and WMAP8 modelled by simply increasing the 
sensitivity by the longer integration time. 



4.1 Dark energy 

Dark energy is, perhaps, the most challenging unsolved problem 
facing modern cosmology. However, as we saw in t]3.4.1l there is 
a strong geometrical degeneracy between 57a, Ok, and w present 
in the CMB. It is therefore necessary to combine CMB information 
with distance measurements in the low redshift Universe in order 
to break this degeneracy and constrain dark energy precisely. 

A promising technique is the measurement of baryon acoustic 
oscillations (BAO) in the galaxy power spectrum. Since the wave- 
length of the BAO is set by the same sound horizon measured in 
the CMB the BAO can be used as a "standard ruler" to measure 
distances. CMB information plays a dual role by calibrating this 
standard ruler and providing a distance measurement at high red- 
shifts. 

Here we consider constraints from a spectroscopic BAO 
galaxy survey of 10,000 sq. deg. spanning redshifts z = 0.5 — 
2. This is a s tage I V space BAO mission in the language of 
Albrecht et al.| j2006) and a possible configuration for the proposed 
JDEM mission. The 68% confidence ellipses in the Q,a — wq plane 
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Table 12. Dark energy constraints for BAO + CMB 

Experiment a(Q_\) cr(wo) relative FoM 



Table 11. Error forecasts for the ACDM + running + tensor spectral index 
model. Second columns shows results for a fiducial model with r = 0.01, 
column 3 refer to r = 0.10. In both cases, we considered the 70 — 143 
GHz channels specifications. 
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Figure 12. 68% confidence ellipses in the SIa — wo plane for BAO in 
combination with CMB experiments. Dashed curve indicates the constraint 
from Planck alone with a weak prior on f! a ■ 



for this experiment are shown in Figure [T2] and summarised in Ta- 
ble [T2j where the constraints have been calculated with a flat prior 
and include representative systematic errors on distance measures 
to each redshift bin and account for non-linear smoothing of the 
acoustic features according to the Fisher matrix pr escription out- 
lined in the Dark Energy Task Force (DETF) report ( Albre cht et al.l 
2006). It is readily seen that the combination of distance measure- 
ments at low and high redshift leads to tight constraints on dark 
energy parameters. The low redshift data is essential in fixing the 
matter content of the Universe (or equivalently h) breaking the ge- 
ometric degeneracy. 

If we expand our parameter set to allow for evolving dark 
energy with equation of s tate w(a) = wo + (1 — a)w a 
dChevallier & Pol arski 2001) the improvement from adding CMB 
can be quantified by the DETF figure of merit (FoM), which is de- 
fined to be proportional to the inverse area contained within the 
Wo — w a error ellipse. These are listed in Table 1 121 WMAP8 pro- 
vides a 16% increase in DETF FoM over WMAP5, while Planck is 
a ~ 2 fold increase in DETF FoM over WMAP5. Thus the goal of 
making precision measurements of dark energy very much require 
precision CMB observations at the Planck level. 
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Figure 13. 68% confidence ellipses in the Qa — f^K plane for BAO in 
combination with CMB experiments. Dashed curve indicates the constraint 
from Planck alone. 



4.2 Curvature 

Since Qk is subject to the same geometric degeneracy as w and 
Qa constraining the curvature is vital for obtaining precision dark 
energy constraints. Typically once constraints at the level of Qk ~ 
2 x 10 are reached this degeneracy is broken dKnox et alj2006t) . 

We apply the same BAO experiment as in the previous sub- 
section to calculate the constraints shown in Figure [T3] The strong 
degeneracy in the CMB data and BAO data have very different 
alignments and in combination lead to constraints at the level of 
J7k ~ few x 10~ 3 . Tablel 1 3 1 summarises the results. 

Including both a variable dark energy equation of state and 
curvature does not degrade these curvature constraints significantly. 
Further, since the curvature is well constrained, we recover con- 
straints on £Ia and w degraded by < 5% from those of Table [T2l 
that were calculated with a flat prior. This illustrates that the degen- 
eracy between curvature and dark energy parameters is broken by 
these experiments. 



4.3 Inflationary parameters 

Constraining inflationary parameters is of paramount importance in 
pointing the way towards a theoretical understanding of first frac- 
tions of a second of the Universe's evolution. The CMB is sensitive 

Table 13. Curvature constraints for BAO + CMB 

Experiment cr(f2A) <j(Qk) 



BAO only 0.05 0.08 

BAO + WMAP5 0.015 0.0037 

BAO + WMAP8 0.014 0.0032 

BAO + Planck 0.0066 0.0017 
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Figure 14. Effect of CMB support on CIP for inflation constraints in the 
n B — n rU n plane. Dashed curve indicates the constraint from Planck alone. 



Table 14. Inflationary constraints for CIP + CMB 

Experiment u(n B ) <r(n run ) relative FoM 



CIP only 
CIP + WMAP5 
CIP +WMAP8 
CIP +Planck 



0.007 
0.005 
0.0049 
0.0026 



0.0032 
0.0022 
0.0022 
0.0015 



1 

1.47 
1.49 
3.07 



( Adshead & Easthej2008l) to quantify this. For this figure of merit, 
both WMAP5 and WMAP8 provide a factor of 1.5 increase over 
CIP alone, while Planck improves the inflationary FoM by a factor 
of 3.0. The combination of Planck+CIP measur es n Iun at the level 
of 10~ 3 (in agreement with lTakada et alj|2006t) , which is roughly 
the largest value consistent with simple slow roll inflation models 
and the current constraints on n s and is also comparable with that of 
some more exotic inflationary models (e.g. Bast ero-Gil et al .120031 : 
IChung et alj|2003h . 



to the shape of the primordial power spectrum over wavenumbers 
ranging from k — 0.002 — 0.2h Mpc -1 , since Silk damping ef- 
fectively erases information about smaller scales. Improving con- 
straints on the tilt n s and running a of the primordial power spec- 
trum will require combining CMB information with information 
derived from the matter power s pectrum, which extends to smaller 
scales increasing the lever arm I Ads head & Easthed feoOS). Preci- 
sion measurements of n s , n run , and r will begin to allow recon- 
struction of the inflaton potential. 

There are currently three techniques proposed for measuring 
small scale power: galaxy surveys, the Lyman-alpha forest, and red- 
shifted 21 cm experiments. The first two techniques are subject to 
systematic uncertainties arising from the non-linear evolution of 
the matter power spectrum, while the latter is still in its infancy and 
may be affected by the imprint of reio nisation on the intergalactic 
medium (IGM) (Furlanetto et al. 2006). Perhaps the most promis- 
ing approach will be high redshift (z — 2 — 6) galaxy surveys, 
which probe the matter power spectrum when the scale of non- 
linearities is small and so may extend to k ~ 2h Mpc -1 or larger. 
For this analysis, we consider the proposed Cosmic Inflation Probe 
(CIP), a space bas ed mission to surve y galaxies in Ha over the 
range z = 2 — 6.5 JMelnick et al. 2004), and assess the importance 
of CMB support for placing constraints in the n s — n Iun plane. We 
assume that CIP will survey 1000 sq. deg. on the sky detecting 200 
million galaxies. CIP is designed to produce an error ellipse for 
inflationary constraints that is comparable to that obtained using 
Planck alone. Since they cover very different scales the combina- 
tion is expected to yield precise constraints on inflation. 

Error ellipses for CIP in combination with CMB experiments 
are shown in Figure [14] Although there is no significant degener- 
acy present in the CMB data, for CIP there is a strong degener- 
acy between n s and n run , since CIP predominantly measures the 
power spectrum on scales smaller than the pivot point fc p i vo t = 
0.05 Mpc" 1 . Planck data provides tight constraints on n s partially 
breaking this degeneracy and leading to much improved constraints 
on n run . 

It is apparent that increasing the observation time from 5 
to 8 years for WMAP shows little improvement in the inflation- 
ary constraints from the combination with CIP, Planck is needed. 
Constraints are listed in Table [14] We may use the generalisa- 
tion of the DETF figure of merit to the inflationary parameters 



4.4 Neutrino masses 

As discussed earlier, the CMB displays a strong degeneracy be- 
tween tJdm and /„ since both affect the energy density in non- 
relativistic matter before decoupling in a similar manner. By adding 
low red-shift information where neutrin os have become non- 
relati vistic this degeneracy may be broken (Lesgourg ues & Pastor! 
2006, for a review of neutrino physics, see). 

In a similar fashion to constraining inflation, constraints on 
neutrino mass benefit from constraints on the matter power spec- 
trum on small scales. In the relatively near future, surveys such as 
the Large Synoptic Survey Telescope (LSST) are likely to sur- 
vey large numbers of galaxy at z < 3 with photometric redshifts. 
These are expected to greatly improve the constraints on neutrino 
masses. A similar analysis taking i nto account the com bination of 
CIP and Planck was performed by |Takada"e t al. (2006). Neutrinos 
show themselves in the matter power spectrum as a suppression 
of power on small scales where relativistic neutrinos are able to 
freestream out of gravitational potentials smoothing the distribution 
of matter. If neutrinos have mass then at late times they become 
non-relativistic and begin to cluster in the gravitational potentials 
set by the dominant dark matter component. Since neutrinos of dif- 
ferent masses would become non-relativistic at different times the 
amplitude of power on small scales as a function of redshift is a 
useful probe of neutrino mass. 

Figure Q3] shows constraints in the Wdm — fv plane for f u = 
0.01, corresponding to a total neutrino mass of M v = 0.1 eV. 
We assume that LSST measures galaxies with photometric red- 
shifts (with uncertainty o z = 0.04) in six redshift bins over the 
range z — 0.5 — 3. We model the observed density of galaxies as 
n g (z) — 640,z 2 e _z / ' 35 arcmin -2 . Note that the uncertainty is to 
some extent driven by the photometric redshifts. If LSST achieved 
a z = 0.02 then the constraints would be significantly improved 
and most of the constraining power would come from LSST alone. 

LSST alone does a reasonable job of measuring f v , but by 
placing a tight constraint on oj& m it also breaks the degeneracy 
present in the CMB data leading to a much tighter constraint on 
fv Once Planck is added to LSST a constraint on M v at the level 
of 0.05 eV is achievable - sufficient for a detection and comparable 



5 http://www.lsst.org/Science/DETF.shtml 
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Figure 15. Effect of CMB support on LSST for neutrino mass constraints 
in the u)dm — fv plane. Dashed curve indicates the constraint from Planck 
alone. 



Table 15. Neutrino mass constraints for LSST + CMB 

Experiment cr(u>d m ) <^{fi/) 



LSST only 0.00097 0.04 

LSST + WMAP5 0.00052 0.0082 

LSST +WMAP8 0.00046 0.0075 

LSST +Planck 0.00029 0.0047 



to the value expe cted from the observed neutrino mass splittings 
dFogli et all2 008). 



5 CONCLUSIONS 

In this paper we performed a MCMC estimate of cosmological con- 
straints expected from the upcoming Planck data. The actual con- 
tribution by Planck to cosmology will ultimately depend on a num- 
ber of factors, especially on systematics control and foreground 
cleaning. In this work we adopted a simple approach to foreground 
cleaning. We supposed that a number of frequencies will be used 
to fully remove foregrounds from the remaining channel, which are 
then used to constrain the cosmological parameters, and discussed 
how such constraints depend on the frequencies used. 

As a general result, we find that a combination of 70, 100, and 
143 GHz channels allows Planck to achieve 85 — 90% of its po- 
tential, corresponding to a factor 3-4 improvement over WMAP 
5 year results for the minimal ACDM model. Adding the remain- 
ing LFI frequencies as well as the 217 GHz channel essentially 
accounts for all of Planck capabilities in terms of cosmological pa- 
rameters. Higher frequencies, while useful for foregrounds control 
and additional science, do not significantly increase the accuracy 
on the main parameters. 

For the cosmology currently preferred by WMAP5 data, 
Planck will provide a cosmic variance limited determination of the 
temperature power spectrum up to £ ~ 1000, while the signal-to- 
noise ratio will be unity around £ ~ 1600 — 1800 depending on the 
channels considered. This high-<? temperature information is essen- 
tial for determination of n s and fib (30-40% improvement wrt con- 
straints obtaining including only multipoles up to £ max ~ 800). 
Multipoles in the range 2000 < £ < 2500 only provide a few 
percent contribution to the accuracy on cosmological parameters, 



however this estimate does not consider gravitational lensing and 
effects like beam systematics which introduce correlation between 
different multipoles. The impact of multipoles in this range over 
the final Planck constraints will strongly depend on the relevance 
of these effects. Future experiments, achieving a cosmic variance 
limited measurement of the TT power spectrum up to £ ~ 2500 
will improve Planck constraints on u>b, n s and 9 by ~ 30 — 40%. 
We find that having a beam slightly bigger than that quoted de- 
grades constraints by ~ 10% on n s and fib. However, adopting a 
slightly wrong value for the FWHM in the data analysis does bias 
the results for several parameters significantly. 

Planck will provide a factor of ~ 9 improvement over WMAP 
in determining the upper limit for the tensor-scalar-ratio. Losing 
the polarization information from the 143 GHz channel would de- 
grade this by a factor of two. Being unable to constrain 73-modes, 
as results of foreground contamination would also degrade the er- 
ror on r by a factor of ~ 2, with a slight dependence on the fidu- 
cial value of the tensor-to-scalar ratio. Planck alone will be unable 
to constrain the tensor spectral index tit- Further, if the full 70- 
143 GHz range is available the degeneracy between n s and n Iun is 
broken allowing constraints on both a factor of five beyond those 
achieved by WMAP. Taken together, Planck will greatly improve 
our constraints on the inflationary parameter space. 

Significant degeneracies remain in the Planck data due to in- 
trinsic degeneracies in CMB. The strong geometric degeneracy in 
the angular diameter distance determination remains, although er- 
rors on parameters like fin and w will improve by a factor 2 com- 
pared to WMAP5 thanks to a better determination of the angular 
scale of the CMB acoustic peaks. Neutrino properties are much 
better constrained (especially N e g for which the current degener- 
acy is totally removed), but a strong degeneracy between f v and 
Wdm remains and considering neutrinos in the analysis increases 
the estimated error on as by a factor ~ 8. 

We have verified that all of these intrinsic degeneracies can be 
broken by the addition of extra low-redshift data sets. Large scale 
structure measurements provide a measurement of fiA breaking the 
geometric degeneracy between fiA, w, and fix present in CMB 
data alone. For a JDEM type mission, Planck improves the con- 
straint on w and fiK by almost a factor of two over using WMAP5 
data. Although Planck does not see a strong degeneracy between n s 
and n lm , galaxy surveys like CIP, which probe smaller scales, do. 
Combining CIP with Planck breaks this degeneracy leading to con- 
straints on the running that are a factor of 2 better than from Planck 
alone. Finally, by measuring the non-relativistic matter content to- 
day, galaxy surveys such as LSST break the degeneracy between 
u)m and f v improving constraints on /„ by an order of magnitude 
to the level of the observed mass splittings. On all of these param- 
eters using Planck offers at least a 50 — 100% improvement over 
WMAP5, illustrating its importance as a support for other data sets. 
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